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Abstract 

The anisotropic s = ^ anti-ferromagnetic Heisenberg chain in the presence of an ex- 
ternal magnetic field is studied by using the standard quantum renormalization group. 
We obtain the critical line of the transition from partially magnetized (PM) phase to 
the saturated ferromagnetic (SFM) phase. The crossover exponent between the PM 
phase and anti-ferromagnetic Ising (AFI) phase is evaluated. Our results show that 
the anisotropy(A) term is relevant and causes crossover. These results indicate that 
the standard RG approach yields fairly good values for the critical points and their ex- 
ponents. The magnetization curve, correlation functions and the ground state energy 
per site are obtained and compared with the known exact results. 
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1 Introduction 



Systems near criticality are usually characterized by fluctuations over many length scales. At 
the critical point itself, fluctuations exist over all scales. At moderate temperatures quantum 
fluctuations are usually suppressed compared with the thermal fluctuations. However if 
the temperature is near zero, the quantum fluctuations especially in the low-lying states 
dominate the thermal ones and strongly influence the critical behaviour of the systems. The 
study of the ground state and its energy is thus of centeral importance for understanding 
the critical behaviour of such systems. 

The technique of Renormalization Group (RG) has been so devised to deal with these 
multi-scale problems[l|, 0, ^, ^. In the momentum space RG which is suitable for studying 
continuous systems one iteratively integrates out the small scale fluctuations and renormal- 
izes the Hamiltonian. In the real space RG, which is usually performed on lattice systems 
with discrete variables (i.e. quantum spin chains), one divides the lattice into blocks which 
are treated as the sites of the new lattice[^. The Hamiltonian is divided into intra-block 
and inter-block parts, the later being exactly diagonalized, and a number of low lying energy 
eigenstates are kept to project the full Hamiltonian onto the new lattice. The accuracy of 
the method is determined by the number of states kept and specially is sensitive to the 
boundary conditions!^, ^ § which is considered for the block Hamiltonian. The detailed 
form of this projection in fact differentiates various versions of the real space RG, ranging 
from the standard RG to the recent Density Matrix RG (DMRG)@. Each of these versions 
have their own advantage and disadvantages. 

The Ising model in a transverse fleld (ITF) and anisotropic XY model in a transverse 
fleld (AXYTF) have been studied in |ll| using both standard RG and DMRG methods. 
There, it has been concluded that DMRG gives accurate results for the ground state energy 
and correlation functions in both models, but the standard RG method where the number 
of states kept are not few can give better results in determining the location of critical 
points and critical exponents. In this direction we have been motivated to study a more 
general model, the anisotropic Heisenberg model in the presence of an external magnetic 
fleld (xxz-|-h) by using the standard RG method to compare its results with the known 
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exact ones. This RG study allows us to have analytic RG equations which gives a better 
understanding of the behaviour of the real space RG method at the critical points. We 
have studied the xxz+h model because of its richness in the phase-diagram where there 
are different critical behaviour. In this study we have succeeded to obtain the critical line 
between the partially magnetized (PM) and saturated ferromagnetic (SFM) phases, to a 



good accuracy, compared with the known results[|T2|. We have also derived the crossover 
between the PM phase (small anisotropy — 1 < A < 1) and the anti- ferromagnetic Ising 
(AFI) phase (large anisotropy A >> 1) and calculated its exponent (0 > 0) which verifies 
the relevance of anisotropy to the crossover phenomena. These results which come out of 
an RSRG by keeping only two states in each block confirm that RSRG is a good candidate 
to study at least the qualitative behaviour of the quantum lattice systems in the quantum 
critical region. 

In this letter we have studied the xxz+h model by RSRG method where the block length 
is three (n^ = 3). In the next section we will introduce the xxz+h model and discuss its 
critical behaviour as derived by other methods. In section 3 we will discuss about different 
types of constructing RG-equations and obtain the analytic RG-equations for this model. 
Using these equations we will describe the phase-diagram of this model. We will discuss 
the critical behaviour of the xxz+h model by RG-equations in section 4. In section 5 we 
will compare some of the results i.e. ground state energy and correlation functions with the 
known exact results. The paper ends with a conclusion. 



2 The Model 



The anisotropic spin | anti-ferromagnetic Heisenberg model, or xxz chain is one of the most 
studied quantum spin systems in statistical mechanics. It is also a classic example of one 



dimensional integrable quantum spin systems |T^, 13, 14]. The xxz chain gives us the first 
example of a critical line with critical exponents varying continuously with the anisotropy. 
The model has also been studied by using the conformal invariance idea[|l^, where the 
critical fluctuations along the critical line are governed by a conformal field theory with 
central charge c=l. The RG study of xxz model has been performed by Rabin [p!^ . At 



A = 0, even in the presence of magnetic field, the model can be mapped by Jordan- Wigner 
transformation to free-fermions [[l7| , which is exactly solvable. This system has also been 
studied within the RG formalism [pJ]] . 

The Hamiltonian in the presence of an external magnetic field is 
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H{J, A, h) = + s\s\^, + As^.^, + K) (1) 

where J > 0, A is the anisotropy parameter, which in the anti-ferromagnetic region is taken 
to be greater than or equal to —1, and /i > is the strength of the external magnetic field. 
The effect of a uniform magnetic field in the phase diagram of the xxz chain is to extend the 
critical phase over a finite region which is partially magnetized and delimited by a critical 



line where the chain becomes saturated ferromagnetic fT^, |15[. Uniform external magnetic 
field does not destroy the exact integrability of the quantum chain but the coupled integral 
equations for the spectral parameter do not have closed analytic solutions. Then the only 



results are numerical or perturbative onesM, [15 
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The Hamiltonian H{J^ A, h) is related by a canonical transformation U = e 
to H{J, —A, h), i.e. UH{J, A, h)U^^ = —H{J, —A, h). This gives a relation between the 
antiferromagntic( J > 0) and ferromagnetic (J < 0) cases. At A = 1 and h = the Hamilto- 
nian exhibits an su(2) symmetry. For A 7^ 1, it exhibits a quantum symmetry stig(2)[|l^. If 
h ^ 0, the only symmetry is U(l). Let us now begin the RG study of the xxz-|-h model. 



3 Renormalization Group Equations 

The implementation of RSRG is based on two important points, the size of blocks and the 
number of states kept in each step of RG. Both of them would have significant effects on the 
RG flow. Here we choose a 3-sites h\ock{nB = 3) for the renormalization process. In this 
case the two lowest energy states of the block Hamiltonian preserve the symmetries of the 
Hamiltonian and lead to a self similar Hamiltonian. Moreover at large A and h the level 
crossing of the ground state in the block occurs at a coupling constant which is exactly its 
critical value (this will be explained latter, see also appendix). Finally taking larger blocks, 
renders an analytic RG-equations, difficult to obtain. 
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After dividing the whole chain into 3-sites blocks, the first step of RSRG is to divide the 
Hamiltonian into two parts, the intra-block Hamiltonian (H^) and the inter-block Hamil- 
tonian (H^^). There are several choices for doing this decomposition. In our prescription 
we choose the decomposition which is sketched in fig.l. In this case the block Hamiltonian 
is 

< = J[sts'2 + 44 + 44 + 44 + A(44 + 44) + h{sl + 4 + 4)], (2) 

where represents the block number i.e. = Y^^=iH^- The block Hamiltonian is 

diagonalized exactly and then the two lowest lying states are kept to span the truncated (or 
effective) Hilbert space. Thus the embedding operator(T) is constructed to be 

T=\a){+\ + \P){-\ (3) 

where \a) and \j3) represent the two low-lying eigenstates of and |+), |— ) are the renamed 
base kets for the effective Hilbert space. 

There is a level crossing at /i = /iq (eq.(§)) in the spectrum of H^, where the ground 
state changes from the = ^ state to the = ^ state {S^ = + + Note that Hq 
comes from the finite size effects of a 3-sites block and reaches the critical value of external 
magnetic field (he) by increasing the block sizes (n^ oo). In that case the ground state 
changes from a PM state (m 7^ 0) to a SFM state (|Tn| = 0.5). Thus in the absence of 
magnetic field, the ground state of the block Hamiltonian is a spin | doublet. As h is turned 
on weakly, we enter in a Zeeman regime and the doublet splits into two states. This is true 
as long as h < Hq. For strong magnetic fields, we are in a regime in which the ground state 
is a singlet with all spins down. This should correspond to h > ho- The results of this 
computation follows : for /i < /iq : 



|a) = &| + +-) +a| + -+) +6| -++), (4) 
\P) = -b\-~+)-a\-+-)-b\ + —), 



where 



3A + V^T8 
^o(A) = , (5) 
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a 



2x + 2 
V6 + 12x2 



2x - 1 
V6 + 12x2 



and 



for h > h 



X 



2(A-1) 
+ A + 3VA2 + 8' 



• 



|a) = -b\ +) - a| - +-) - 6| + ), 

m = I---)- 



(6) 



Having the embedding operator at hand, the operators (observables) are renormahzed as 

Qf ^ T^OT. (7) 

By using the above equation one can obtain the renormahzation of operators. Thus for 
h < ho we obtain the following relations : {h < Hq) 



Tts^(3^T = -2abs"' 
Thl^^-^T = 2abs'y 



TUIT = -2b'^s"' 
TU%T = 2bh'y 



(8) 



1(3) 



2n Jz 



2a^)s 



We find the same renormahzation for si and S3 because of the symmetry in site 1 and 3 in 
the block (1 — 2 — 3). We will obtain similar relations for h > Hq : 



-bs'- 



J Si(3)J - I + S 



Th^T = -as"" 
.2 1 



(9) 



a 



In the above equations 3'°" is the effective operator in the effective Hilbert space of the block 
(new sites in the renormahzed chain). By considering the interaction between blocks and 
using the above equations (§),(|^) we will obtain the renormahzation of coupling constants 
in the Hamiltonian : 

J' = Aa%^J, 

\h\ < ho{A) : 



A' = -SiA 



(10) 



h' 



4a262"-5 



J' = 

\h\ > ho{A) : I A'^b'^A, (11) 

_ h' = ^{\h\-ho + ^A)sgn{h). 
The above RG-equations show that the renormahzed Hamiltonian is of the same form as 
the original one. The critical behaviour which can be obtained from these equations will be 
discussed in the next section. 

Let us consider an extreme case where both h and A go to infinity and J goes to zero 
such that JA, Jh and remain finite. In this case the Hamiltonian reduces to a simple 
anti-ferromagnetic Ising (AFI) model in the presence of an external magnetic field, which 
shows a first order transition from a classical anti-ferromagnetic(Neel) ordered phase (m = 
0, sm — jj: l)'sf = 0.5) to a saturated ferromagnetic phase (SFM) (|m| = 0.5, sm — 0). 
We can write this Hamiltonian as 

HAFi^kY^islsUi + gst) (12) 

i 

where 

k^JA>0 , g = ^ (13) 

At large A (A >> 1) we have /iq — A ,a ^ 1 and b ^ . Then the RG-equations 
reduces to the following equations. 



\h\<A : 



46V, 



J' 
A' 

~ 462"'' 



(14) 



\h\> A : 



(15) 



J' = b^J, 
A' = b^A, 
^ h'^^{\h\-A)sgn{h). 
These RG-equations give exactly the critical point gc — ^ and the ground state energy of 
the AFI model(see appendix), which will be discussed in the next section. 
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4 Critical Behaviour 



In this section we analyze the RG-equations and its critical behaviour. The phase diagram of 
the obtained RG-flow(eqs. (p!0| , pTD ) is depicted in fig. 2. This phase diagram consists of three 
different phases, partially magnetized (PM), classical anti-ferromagnetic (AFI) and saturated 
ferromagnetic (SFM) phases. 

There are five fixed points in the phase diagram, 
(i) XX represents a spin ^ xx model without external field, 

(m) XXTF is the critical point of the xx model in the presence of a transverse field, 

[iii) lAFH represents the critical point of the xxz model in the absence of an external field, 

(iv) AFI represents a classical anti-ferromagnetic Ising model with a long-range Neel order, 
(w) SFM represents a saturated ferromagnetic phase where all spins align in the direction of 
the external field. 

In the SFM phase the RG-fiow has a well defined behaviour and goes to the SFM fixed 
point for any value of h > hc{A)- But when h < hc{A) and — 1 < A < 1 the RG-fiow 
represent a massless phase in which J*^") and A*^") — > in the limit n oo {n is the 
number of RG steps). The RG fiow in the PM phase has a cyclic nature. Since it refiects 
a sequence of level crossings between states with different values of the total 5"^ induced by 
varying the magnetic field. The recurrence of this level crossing in the process of RG leads 
to the oscillatory behaviour of the RG fiow. If we imagine of a 3-dimensional RG fiow, its 
projection onto the h — A plane will have some closed paths. However if the initial point is 
in the PM(or SFM) phase, it will go to the XX(or SFM) fixed point finally. Therefore we 
conclude that the RG fiow in the PM phase can be sketched as in fig.2. At the end of this 
region on the A = line there exists a fixed point at he = 0.943 which separates the PM and 
SFM phases. The eigenvalues of RG-fiow at this critical point are given in table-1, which 
give a relevant direction on the A = line and an irrelevant direction along the critical line 
he = 0.943A -|- 0.943. This critical line is obtained from the behaviour of the correlation 
functions. We have calculated (0|s^s|'_|_]^|0) and plotted it versus h (fig. 5) which shows the 
entrance to the SFM phase for h > hc{A) in the — 1 < A < 1 region. It is an interesting 
result which is obtained by a 3-sites block RG and can be compared with the exact result 
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/ic = A + 1 Although the obtained critical value for hc{A) has a slight difference in the 
coefficients but preserves the linear form of the critical line. Our data for the ground state 
energy (eo) show that hc{A) represents a critical line of 2nd order transition in which Co and 
%^ are continuous at /ic(A) and is confirmed by the analytic results for eo{h) at A = [Tl\. 



Our results along the A = and h = lines recover the results of Drzewinski ||Tl[] and 
Rabin |]16| respectively. 



For the XXTF fixed point we have calculated the critical exponents which have been 
written in table-1. If R{h) represents the renormalization of h along the A = line the 
correlation length exponent (u) is given hj u = i^^^^T^flyy- The dynamical exponent (z) is 

ln{ — ) * 

z = ■ The critical exponent a connected with the specific heat is calculated from the 

hyper-scaling relation 2 — a = d*v [|T^, where d* = d + z {d is the spatial dimension). The 
critical exponent {3, related to the magnetization is given by /? = in(ji^h*) - These results 
show good agreement with the exact ones. 

The other phase in the phase diagram is a classical anti-ferromagnetic phase. Let us first 
look at the exact solution of this model with the Hamiltonian as in eg . {\L2^ . By a simple 
argument we can find that the ground state is a Neel ordered state whose energy per site is 
^ = Co = ^ for < 5f < 1 and is a saturated ferromagnetic state for > 1 where Cq = 
(see appendix). These values for the ground state energy(eo) show a discontinuity of ^ at 
(7=1, which means the phase transition at this point is classified as a 1st order transition. 
By using the previous definitions for g = the RG-equations in (|T^ and (^) give a fixed 



point at (7 = 1 in the limit A ^ oo, which is equal to the exact critical point. The RG-fiow 
in (0) and (|15|) show a fixed line at A — > oo for all g < 1. This means that there is a unique 
ground state for any value of < (7 < 1 and the distinction between two different g values is 
only due to the excited states of the Hamiltonian. 

We have also calculated the ground state energy by using the RG-equations (|T4[) and 
(P!3|) in a hierarchical way by accumulating the energies of the blocks (see appendix). The 
obtained result is equal to the exact result for the ground state energy. 



(16) 
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Thus the hmiting case of our RG-equations at h,A >> 1 exactly describe the classical 
anti-ferromagnetic Ising model. 

There is an interesting point in the phase-diagram. When — 1 < A < 1 and h < hc{A), 
the model represent a PM phase with —0.5 < m < and undergoes a 2nd order transition 
to the SFM phase(m = —0.5) at /ic(A)(m is continuous at the transition point). But at 
/i, A oo; < 1 when the model represent a Neel ordered phase (AFI) with m = and 
sm = 0.5, it undergoes a 1st order transition to the SFM phase(m = —0.5) at /i^ = A (m 
is discontinuous at the transition point). We have examined the continuity of the ground 
state energy and its derivatives up to the third order numerically at A = 1; /i 7^ 0. Therefore 
there is a continuous change in the critical exponents by increasing A when h < he- This 
shows a crossover between PM and AFI phases. This change in the universality class is due 
to the reduction in the number of components of the spin operator (i.e. three component 
3^,3^, 3^ in the PM phase and effectively one component 3^ in the AFI phase). The crossover 
exponent (p = ^ has been calculated to be 0.63 which verifies that the coupling A is relevant 
and causes crossover (Aj = n^^ ; where Aj is the eigenvalue at lAFH fixed point). 

One can also retrieve the critical line 0<A<1; /i = Oby RSRG using the 3Uq{2) 
symmetry [pn| of the xxz chain which differs only at the ends and is not important in the 
thermodynamic limit (A^ —>■ 00). However it is not suitable for A > 1 case and also in the 
presence of the external field {h) this RG prescription could not describe the critical line 

5 Energy and Correlation Functions 

In this section we describe some more results which have been obtained by RG-equations in 
section-3. In fig.3-a we have plotted the ground state energy per site(eo) of the XX(A = 0) 
model versus external magnetic field(/i). We have also compared the RSRG results with the 
exact one ||1 7| . In the PM phase (fig.3-a) there is a discontinuity in Cq which is due to level 
crossing (finite size effects) of 3-sites block. This level crossing occurs as ho{A = 0) passes 
the value 0.707. Thus we have not considered this point as a critical point. Fig.3-a shows 
good agreement with the exact results in the SFM phase(/i > h^) with a slight difference in 
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the PM phase. It has been shown[|, ^ ^ that this difference is due to boundary conditions 
in an isolated bloclc which neglects the remaining part of the chain. We have introduced 
a modified scheme to decrease this difference in the h = case [§. By using the RG- 
equations of section-3 the ground state energy for different values of A can be calculated 
which will show the same behaviour as in fig.3. We have plotted the ground state energy for 
A = — 1, —0.5, 0.5, 1 cases in fig.3-b. The ground state energy at A = —1 in fig.3-b confirms 
that the model represents the ferromagnetic Heisenberg model in the presence of an external 
magnetic field, where its ground state energy is proportional to the strength of the external 
field. 

We have plotted in fig.4-a, the magnetization(m) versus external magnetic field(^) for 
A = 0. It has been compared with the known exact result, which shows good agreement 
qualitatively. The step form of the RSRG results in this figure is due to the cyclic nature of 
the RG-equations in the PM phase. This is related to the nature of the anti-ferromagnetic 
problem. The magnetization curve reflects a continuous sequence of level-crossings between 
states with different values of the total induced by varying the magnetic field. The 
variational ground state which is obtained here is however owing to the highly degenerate 
energy level structure. The recurrence of this level crossing in the process of RG leads to 
the oscillatory behaviour of the RG flow. This oscillation is trapped by a metastable state 
which leads to a jump in the magnetization curve. Fig.4-a confirms that for h < he the 
model is PM where m 7^ and reaches the SFM phase (m = —0.5) at h = h^. We have 
also plotted the magnetization versus external field in fig.4-b for A = —0.5, 0.5 and fig.4-c 
for A = 1.1, 1.5 which show a similar behaviour as in the A = case, but the critical point 
where the model becomes SFM is different. 

One of the important quantities which can be used to show the critical behaviour is the 
z-component of spin-spin correlation function. We have plotted (0|s^s^_,_^|0) versus external 
magnetic field(/i) in fig. 5 for different values of anisotropy parameter A. Transition from the 
PM phase to the SFM phase occurs at two steps. The first jump is due to level crossing at 
hQ{A) which is not the critical point and does not terminate at the SFM phase. But the 
second jump which achieves the SFM phase ((0|sfs^'_,_]^|0) = 0.25), corresponds to the critical 
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point of the transition from the PM to the SFM phases. The location of critical point is at 
the same point in which the magnetization curves (fig.4-a, fig.4-b and fig.4-c) reaches the 
saturated value (m = —0.5). The critical point hc{A) for different values of A which has 
been obtained in fig. 5 yields the linear relation hc{A) = 0.943A + 0.943. This result shows 
that RSRG is a good candidate to describe the critical behaviour of quantum systems |1TT||. 
We have also calculated the magnetization(m) in the middle range A > 1 and found that the 
critical line hc{A) = 0.943A + 0.943 is also valid in this region. Since the RSRG method is an 
approximate scheme there is a small discrepancy between the obtained critical line and the 
exact one {he = A + 1). This error is related to the small isolated blocks which are considered 
in this method. In other words the quantum fluctuations of a highly correlated system in 
a large lattice can not be simulated by a small number of eigenkets of an isolated block. 
However as A ^ oo and the model becomes classical one, exact results can be obtained by 
RSRG method (see appendix). This describes the discrepancy in the limit A ^ oo of hc{A) 
which is (7c = ^ = 0.943 + 0{^) and gc = ^ which can be obtained by RSRG method in the 
limit A — > oo of the initial model. We have also calculated the z-component of spin-spin 
correlation function in terms of distance. In fig.6 we have plotted (0|s^s^^_^|0) at A = 1 
for different values of h below and above its critical point (/ic = 1.886). When h = 0, the 
correlation length is small(in the order of lattice spacing) and the correlation function goes 
rapidly to zero after a few lattice spacing. As h increases to its critical value the correlation 
function becomes nonzero for long distances and shows exactly the SFM phase above the 
critical point (/;, > he)- The same behaviour has also been observed for other values of the 
anisotropy parameter A. 



6 Conclusion 

We have considered the anisotropic anti-ferromagnetic Heisenberg chain in the presence of 
an external magnetic field by RSRG. We have sketched the phase diagram of this model 
for A > — 1 and h > in fig.2, the phase diagram for A > —1 and /i < is the mirror 
image of the previous case. We have obtained three distinct phases in the phase diagram. 
The partially magnetized (PM) phase with m 7^ 0, the saturated ferromagnetic (SFM) 
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phase with m — —0.5 and the Neel ordered phase (AFI) where m — 0. By computing the 
magnetization and the z-component of spin-spin correlation function we have calculated the 
critical line hc{A) = 0.943A + 0.943 between the PM and SFM phases which causes a 2nd 
order transition. But at A — > oo the transition from AFI to SFM phases is 1st order at 
the critical point gc — i- We have observed that for h < he increasing A causes a crossover 
between PM and AFI phases which changes the universality class of the model. The crossover 
exponent has been calculated to be = 0.63 which confirms the relevance of the anisotropy 
A in the crossover phenomena. 

By using the analytical RG-equations we have obtained the critical exponents at the 
XXTF fixed point. Although the obtained critical exponents are not accurate compared 
with the exact results they show good agreement with them. The ground state energy and 
correlation functions calculated in the PM phase, show qualitatively good results, but some 
discrepancy due to the boundary conditions of the isolated block in the RG procedure is 
present. However all the results in the SFM phase are completely accurate, because the 
ground state of the whole chain in this phase is a simple juxtaposition of the ground state of 
the isolated blocks and there is not any boundary condition effects for an isolated block as 
in the PM phase. In the AFI fixed point when the model is a classical one, the limiting form 
of our analytical RG-equations (A — > cx)) give the exact results for the ground state energy 
and the critical point gc — 1- We have shown that the critical line /ic(A) is also valid in 
the middle range A > 1. Finally we conclude that the standard quantum RG (RSRG) gives 
qualitatively good results for the critical behaviour of the system. However the quantitative 
results for the location of the critical point and critical exponents is much better than its 
results for ground state energy and correlation functions with respect to the known exact 
results. 
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8 Appendix : Classical anti-ferromagnetic Ising model 
in an external field 

8.1 Exact ground state 

The classical anti-ferromagnetic Ising model in the presence of an external magnetic field is 



given by equation (|T^ 



N 

Hafi = kY^^s^sl^^ + gsl) 
1=1 

where k = J A > and g = ^ is the strength of external magnetic field. We assume periodic 
boundary condition s|^r^_^ = sf. Let us write the Hamiltonian in terms of Pauli matrices, 

u N N 

Hafi = j\Ei^!^Ui " 1) + '^9Y.^^ + (17) 

^ i=l 1=1 

The energy (E) of this chain is 

E=^[-2nj + 2gns + N], (18) 

where the number of up spins, rim is the number of down spins {np + Um = 

N) and rij is the number of boundary walls of flipped spins. The maximum value of is 
obtained by a Neel ordered state {up = rim = Y^^f ~ ^)' ^'^^ arbitrary value of 
it can be written as 

{nf)max = N - \ns\. (19) 
By using the definition of rig, {nf)max can be written in the following form 

2n n < — 

{nf)max = N-\N-2nm\ = { " (20) 

2{N-nm) nm>f 

The minimum value of E will be obtained if ris has its minimum value and rif has its 
maximum value. 

Eo = rmn{E) = ^[N + 2g{N - 2r2„) - 2{nf)max]. (21) 
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When Um < we have j§ = = ^ which is greater than all the energies in the Um > y 
case. Therefore we will investigate Eo in the Um > y region by minimizing E with respect 
to Um, which is 



1-23 
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(7>1 (n„ = A^) 



It is obvious from (p4 ) that for any value of (7 < 1 the ground state energy is due to a Neel 
ordered state (n^ = ^^rij = N) and for (7 > 1 the ground state is a saturated ferromagnetic 
state {rim = N,nf = 0). 

8.2 Renormalization group equations 

At large A(A ^ 00), the renormalization group equations (|T^ and (|T3p can be written for 
k = J A > and g = ^ > in the following form 

f k' = k, 

g<l : (23) 

[ 9' = 9, 

k' = b^k, 

g>l : ' (24) 



The RG equations in (|23| ) give no running of coupling constants and lead to a fixed line 



< g < 1. As far as the ground state energy is concerned, there is no distinction between 
any arbitrary value of < g < 1. But when g > 1 the only fixed point is 5'* = 1 which is at 
the end of the fixed line < g < 1. Thus g* = 1 is the critical value of g which separates the 
g < I and 9 > 1 phases. To be more rigorous and define these phases we will calculate the 
magnetization (m) and staggered magnetization (sm) in these regions. Equations (^ and 
(P) in the limit A^oo {a l,b ^ 0) will be written as 



g<l ; Thl^,^T = s'' ; T^s^^T = -s'^ (25) 

-/ 
2 



9>l ; Ttsf(3)r = — ; Tts^T 



The magnetization (m = jrY^i^ (Oki|0)) is 



m = j^T.m4, + 4, + sim (26) 



N/3 
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1 m 

where |0') is the ground state in the renormahzed Hilbert space (T|0') = |0)). By using 
equation (|25|), the magnetization is calculated to be 

m = -m! ; fi' < 1 (27) 

-11/ 
m = — + -m ; g>l 

where m' is the magnetization in the renormalized chain with y sites. However in the 
thermodynamic limit (A^ oo), m and m! will be equal. Thus equation (^) gives 

m = ; g <l (28) 
m = —0.5 ; > 1- 

Similarly, we can use the definition of staggered magnetization (sm = j^Ylf (0|(— l)'sf |0)) 
and repeat the steps in calculating m, the staggered magnetization is obtained to be 

sm = 0.5 ; g <1 (29) 
sm = ; (? > 1- 

These results confirm the Neel ordered state in the g < 1 region and the saturated ferromag- 
netic state in the g > 1. Note that these results are the same as the exact ones which have 
been obtained in the last section. Had we taken even size blocks for the renormalization 
procedure, we could not obtain these values. 

The calculation of the ground state energy is easily done by accumulating the energy of 
blocks in a hierarchical way. The renormalized Hamiltonian is 

HN{k ,g ) = THN{k,g)T = i 30 

[ '-^^ + k{g-l)j:^/'sr ■ g>l 

Therefore the ground state energy per site is calculated to be 

^ ^ i(f)(l + | + | + -) = X 9<l 

^° kN I 1^ _ 2(g-l) ^ l-2g . 

I, 12 12 4 1 y ^ 

which is equal to the exact ground state energy. This result can not be either obtained by 
taking an even size block in the RG procedure. 



(31) 
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Tables 



Table- 1. Eigenvalues and critical exponents at the XXTF fixed point, both RSRG and 
exact results. 





Ai 


A2 


P 


V 


z 


a 


RSRG 


0.250 


4.000 


0.792 


0.792 


1.262 


0.208 


Exact 






0.5 


0.5 


2 


0.5 
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Figure Captions : 

Fig.l) Decomposition of lattice into block and inter-block part, and different types of 
intra-block (H^) and inter-block (H^^) interactions. 

Fig. 2) Phase diagram of the xxz model in the presence of an external magnetic field 
{h). Filled circles are the fixed points and arrows show the direction of flow. The solid line 
which passes through the (A = —1, h — 0) and(A — l^h — 1.886) points is the critical hne 
he — 0.943A + 0.943. The dotted hne for A > 1 show qualitatively the crossover region. The 
double solid hne at A = oo is the fixed hne < g < I. 

Fig.3-a) Ground state energy per site versus external field (/i) for A = , both RSRG 
and exact results. 

Fig.3-b) Ground state energy per site versus external field (h) for A = —1, —0.5, 0.5, 1, 
which has been obtained by RSRG. 

Fig.4-a) Magnetization (m = -^JliLi (Okl |0)) versus external field (h) for A=0, both 
RSRG and exact results. 

Fig.4-b) Magnetization {m — J2iLi (0|s||0)) versus external field (h) for A = —0.5, 0.5, 
which has been obtained by RSRG. 

Fig.4-c) Magnetization (m — jj: Y.^i |0)) versus external field (/i) for A = 1.1, 1.5, 
which has been obtained by RSRG. 

Fig. 5) z-component of spin-spin correlation function for different value of anisotropy 
parameter (A) versus external field {h). 

Fig. 6) z-component of spin-spin correlation function at A = 1 versus distance (r) for 
different value of external field (h). 
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